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Abstract 



We perform the complete bosonization of 2+1 dimensional QED with one 
fermionic flavor in the Hamiltonian formalism. The fermion operators are explic- 
itly constructed in terms of the vector potential and the electric field. We carefully 
specify the regularization procedure involved in the definition of these operators, 
and calculate the fermionic bilinears and the energy - momentum tensor. The al- 
gebra of bilinears exhibits the Schwinger terms which also appear in perturbation 
theory. The bosonic Hamiltonian density is a local polynomial function of Ai and 
Ei, and we check explicitly the Lorentz invariance of the resulting bosonic the- 
ory. Our construction is conceptually very similar to Mandelstam's construction 
in 1+1 dimensions, and is dissimilar from the recent bosonization attempts in 2+1 
dimensions which hinge crucially on the existence of a Chern - Simons term. 
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It is well known that, in 1+1 dimensions, fermion field operators can be constructed in 
terms of local bosonic fields. This bosonization procedure is of great theoretical interest, 
since it provides a method of mapping different local quantum field theories with a 
priori different Hilbert spaces onto each other. This mapping has proved to be very 
helpful in the analysis of different 1+1 dimensional models. After the pioneering works 
of Coleman |lj and Mandelstam ||, bosonization (and its nonabelian generalization ||) 
has been used in a variety of contexts: relativistic quantum field theory, condensed 
matter systems, string theory. 

Clearly, the extension of bosonization to higher dimensions is highly desirable. There 
were several early attempts f|], but those did not lead to local bosonic theories. Until 
recently it has been tacitly assumed that such an extension is impossible. The feeling 
was that 1+1 dimensions is a very special case since there is no spin in 1+1 dimensions, 
and therefore no "real" difference between bosons and fermions. 

The renewed interest in this problem was triggered by a possible connection between 
the phenomenon of Fermi - Bose transmutation in 2+1 dimensions and novel condensed 
matter systems, most notably high T c superconductors. Polyakov argued [|] that the 
addition of a Chern - Simons term || to the Lagrangian of QED 3 changes the statis- 
tics of charged excitations. His argument was elaborated further in many works j7|, 
and the explicit realization of this Chern - Simons induced mechanism of Fermi - Bose 
transmutation in lattice theories has been given in ||, ||. 

There is, however, an important conceptual difference between the Mandelstam con- 
struction in 1+1 dimensions and the Chern - Simons bosonization in 2+ldimensions. In 
the former case bosonization is achieved essentially through a duality transformation: 
one constructs fermionic operators ip, which carry a global U(l) fermion number charge, 
in terms of a local bosonic field 0, which is itself neutral. The U(l) charge, when ex- 
pressed in terms of 0, is a topological, rather than a Noether, charge. In the 2+1 - 
dimensional construction described above the fermionic field ip is constructed in terms of 
a bose field 0, which carries the same global quantum numbers. The element of duality 
is notably missing here. Moreover, in this picture, a free fermion is bosonized into a 
particle interacting with a vector potential with the Chern - Simons action. This bose 
field is not gauge invariant, and after gauge fixing becomes nonlocal. So one does not 
really construct fermionic operators in terms of local bosonic fields. 

The 2+1 dimensional construction depends crucially on the existence of the Chern 
- Simons term. The Fermi - Bose transmutation occurs only for a fixed value of its 
coefficient. For other values of this coefficient, the procedure gives nonlocal anyonic 
fields with fractional statistics. Another unsatisfactory feature of the construction is 
that it has been consistently implemented only for nonrelativistic fermions [0, [pj. When 



applied to continuum relativistically invariant theories |1(| it suffers from regularization 
ambiguities, and it is not clear how to interpret the formal results. In fact, the covariant 
Dirac field has not yet been constructed, even formally, in this framework. 

In this letter we take a different approach to bosonization. It is the direct extension of 
Mandelstam's procedure to higher dimensions. Our aim is to construct a Dirac doublet 



2 



of fermionic operators in 2+1 dimensions in terms of local bosonic fields. As a working 
hypothesis we assume that, in terms of the bose fields, the fermion number charge is 
topological, and thus the associated current is trivially conserved. We are not aware of 
any theorem that states this, but such is the case in all of the known examples where 
either exact bosonization can be performed (in 1+1 dimensions) jl], 0, or fermionic 
excitations are known to exist in the spectrum of a local bosonic theory in 2+1 [11], or 
3+1 [TT2] dimensions. 



Given this assumption, the theory that suggests itself as a convenient candidate for 
bosonization is QED with one two - component Dirac fermion. The reason is that 
Maxwell's equations, 

r = -d^ (i) 

tell us that, if one takes as basic variables the electric field Ei and the vector potential 
Ai, fermion number is trivially conserved, due to the antisymmetry of F^. The fermion 
number charge is topological, since it is equal to the surface integral of Ei at spatial 
infinity. Since there are no additional flavor symmetries here, it is the simplest model of 
its kind. 

Our strategy is as follows. Consider the Hamiltonian of the theory in the temporal 
(Weyl) gauge A = 0, 

H = ^E 2 + ^B 2 + ifj-f i {id i + eA i )tfj + 7mljtfj (2) 
together with Gauss's constraint, 

diEi = eifj^ (3) 

We will solve the constraint by constructing the doublet of anticommuting fermionic 
operators ip a in terms of the local fields Ei and their canonical conjugate momenta Ai, 
[E l (x) , A* (y)] = i5^5 2 (x — y). Substituting those back into the Hamiltonian, eq.(0), 
we will obtain the completely bosonized form of QED 3 , defined on the physical Hilbert 
space. We will then calculate the fermionic bilinears, and check, that they satisfy the 
tree level algebra, modified by Schwinger terms. The appearance of these Schwinger 
terms is also seen in perturbation theory at the one loop level. 

To construct the operator if> we must first fix the gauge freedom associated with the 
time independent gauge transformations, generated by the Gauss's constraint. We do 
this by considering ip in the Coulomb gauge. When written in terms of these operators, 
the covariant derivative in eq.(^) contains only the transverse part of the vector potential 
Ai. Furthermore, we will work with the gauge - invariant operators, 

^ G {x) = ^ a (x)e ie I d2yeiiy - x)Aiiy) (4) 

where, 

eAx) = -4 (5) 
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is the electric field of a point charge. (We will drop the superscript CG from now on). 
In addition to solving the constraint, eq.@, the fermionic operators ip a must satisfy the 
following conditions^: 

i. Carry unit electric charge, [ip a (x),d i E i (y)] = eip a (x)5 2 (x — y); 

ii. Transform correctly under rotations. For convenience we choose the basis of Dirac 
matrices in which the rotation generator is diagonal, 7 = a 3 , 7 1 = icr 2 , 7 2 = —ia l . 
Then, 

Vi-e^Vi; V2-e-**/V2 (6) 
where is the rotation angle; 

iii. Fermionic bilinears must be local operators, i.e., [ip\ L ipp{x),0{y)\ = for x 7^ y 
and for any local gauge invariant operator 0(x). 

To satisfy the first condition, we take the following ansatz, 

i) a {x) = kAV a (x)e ie f d2ye ' {y - x)Mv) U a (x) (7) 

where A is the ultraviolet cutoff, k is a constant, and the operators V a and U a commute 
with the charge density creation operator. The following form of these operators (which 
is not final!) ensures the anticommutativity of the ip's, and also has correct rotational 
properties, 

Vl ^ x ) = _ w ± J dh(e( x -y)-^E l{ y) . v 2 (x) = -iVl(x) (8) 
U^x) = e -£Jd 2 y0(y-x)diEM . = jjl^ 

Here 9(x) is the polar angle of the point x (for the precise definition of 6{x) see [pp. 
With these definitions we find, 

M*)Mv) = My)Mx)e- my - x) - e{x - y)) = -My)M%) (9) 

= 4(y)M^> mv ~ x) - eix ~ y)) = -4(v)Mx) 
M*)Mv) = My)Mx)e my - x) - e{x ~ y)) = -Mv)Mx) 
M*)4(y) = 4(y)Mx)e- my - x) - e{x - y)) = -4(y)M x ) 
M x )4{y) = 4(y)M%)e~ lw = -4{y)M x ) 

ipi(x)ip 2 (y) = ip 2 (y)ipi(x)e™ = -iki(y)1>i(x) 

The factor kA in eq.(0) takes care of the correct dimensionality of the fermionic fields, 
and ensures the correct normalization of the anticommutators, 

{4(x), My)} = i4(x), My)} = s 2 {x - y ) (10) 

However, the above manipulations are formal, since the commutators of V and U 
with the A - dependent factor in eq.(0) are singular. The expression for the fermionic 

1 Bosonization of the massive Schwinger model in 1+1 dimensions can be formulated in precisely the 
same terms. 



4 



operators needs therefore regularization. We shall regularize it by point splitting. We also 
modify the expressions for V(x) and U(x), so that they become creation and annihilation 



operators of a magnetic vortex of half - integer strength |T3 



V±(x) = —iexp 
Ui(x] 



2^ 



d 2 y \(9(x -y)- rfdtEiiy) + 2TxG^ 2 \y - x)e lJ d i E j {y) 



exp< - 



2e 



d 2 y [6{y - x)d i E i (y) + 2^ 2 \y - x)e ij d i E j (y) 

V 2 (x) = -iV}{x); U 2 (x) = U\{x) 

Here G^ 2 \x — y) = — 4- ln(x 2 /i 2 ) is the two - dimensional Green's function of the free 
massless field with IR cutoff \i. The physical interpretation of the operators V a and U a 
is clear from the following commutation relations, 



[V 1 (x),B(y)} = -* Vl (x)6 2 (x-y); [U l {x),B{y)\ = -U l {x)8 2 {x - y) (12) 



The regularization we employ does not change either the anticommutation relations or 
the rotational properties of the operators, but is crucial to ensure the locality of the 
bosonized forms of the fermionic bilinears. 

The regularized expression for the fermionic operators involves an average over the 
direction of the regulator rj, which implements the point splitting, 



kA 



kA 



i)i(x) = lim — / dfje 1 2 ipi(x); ip2{x) = lim — / df]e % 2 i>%{x 



A^oo 2TX 



where 



^(x) = V 1 {x + ^I^^-^^U^x - 77); 
$Z(x) = V 2 {x + v )e ie f d2ve * {v - x)A > iy) U 2 {x - rj) 



(13) 
(14) 



The phase factors introduced into the averages ensure the correct rotational properties 
of the fermionic operators. The length of the regulator rj is taken to be proportional to 
the inverse of the UV cutoff, \rj\ oc 1/A, the integral is over the angle of the vector rj, 
and 17 is the unit vector in the r] direction. It can be explicitly checked that, 



{tp a (x),ijp(y)} = 0, \x-y\»l/A 



(15) 



Therefore, in the limit A — > 00, we regain the standard anticommutation relations. 
Equations (|ITD , flT3| ) and ( H[) are our final expressions for the fermionic operators in 
terms of the bose fields E i and A4. 

We next calculate the fermionic bilinears. Those should also be defined using the 
point splitting procedure. The following definition of the bilinears is a natural one, 



Jr(x) = ip(x)Tip(x) 



57T 



dee ixr(i) I \ip\x + e), 7°r^(x - e) 



Jl-f 1 > 



|e|oc|rj| 



(16) 
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It is implicitly understood in eq.(|16D that the limit A — > oo is taken after (independent) 
averaging over the directions of e and 77. The insertion of the Wilson factor is appropriate 
for the definition of bilinears in a gauge theory. Since we are constructing the Coulomb 
gauge fermions it is only the transverse part of the vector potential, A tr , that appears in 
the Wilson factor. The phase exp{ixr} is inserted to project onto the relevant irreducible 
representation of the 2D rotation group while averaging over e. Thus for T = 70 and T = 1 
we have \( e ) — 0, while for r = 7+ = 71 + i^2, x( e ) — 0(e) an d f° r T = 7- = 7i — ij2, 
x(e) = —0(e). The ratio r]/e is arbitrary, but the final results should not depend on it. 
We also normal - order the bilinears with respect to the perturbative vacuum. 

We have calculated the bilinears using the definitions of the fermionic operators given 
above. The calculation is performed, in analogy to the 1+1 dimensional case, by expand- 
ing all quantities in powers of the inverse cutoff. This procedure has a certain caveat 
that one has to keep in mind. The operators which are multiplied by inverse powers of 
A, and are formally small, may, in fact, give finite contribution in the continuum limit, 
if these operators have high enough dimensions. This problem arises also in 1+1 dimen- 
sions, where the effect of the higher order terms is to renormalize the coefficients of the 
lower dimensional operators in a way consistent with the symmetries of the problem []. 
Furthermore, in order to perform the actual calculations, we must control the ultraviolet 
behavior of the bosonized theory. This behavior is different from the naive one, obtained 
in perturbation theory. This is related to the fact that, apart from the ultraviolet cutoff 
A, the bosonized theory has an additional ultraviolet scale /i = (e 2 A 2 ) 1//3 . The appear- 
ance of this scale can be seen as follows. The fermionic operator defined in eq.([13|) is 
essentially a product of a vortex and an antivortex operator, separated by A -1 . The 
vortex operator scales, at short distances, as an exponential, 



< V(x)V*(y) >oc \x 





' c 


\x — y \ a exp 


~^\ x - y\ 



1 2-a 



(17) 



where c is a constant, and a/2 is the scaling dimension of the electric field, < E(x)E(y) >oc 
7. Assuming the perturbative scaling of the electric field at short distances, a = 3, 



\x-y\ 



we find that the fermion operator that we have constructed, at distances larger than 
1/A, but still small relative to any physical distance scale, behaves, roughly, as, 



< tjr (x)ij}^(y) >oc \x — y^exp 



cffc 



e 2 A 2 \x — y\ 3 

The scale /1 is therefore clearly a crossover scale. At distances larger than, 



\x — y\ oc 



e 2 A 2 



fJL" 



(19) 



the exponential in eq.flTj) can be expanded in power series. In that case the fermion 
propagator scales as a power, which is consistent with its perturbative behavior. 



2 As in 1 + 1 dimensions, this renormalization ambiguity affects only the spatial components of the 
current. 
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At distances smaller than the non - point - likeness of if) becomes important. 
As a result the scaling behavior at these short distances in the bosonized theory must 
be different from the perturbative one. In fact, if the fermionic operator is still to scale 
with a power law for /j, << 1/x << A, the leading UV behavior of the propagator of 
the electric field must be < Ei(x)Ei(y) >\ x -y\<i/fi~ \x- y \ 2 • Dimensional considerations, 
utilizing the bosonized Hamiltonian 113 are consistent with this result. 



We thus assume the following form of the UV asymptotics of the correlators of the 
bosonic fields, which is the most general one consistent with this UV scaling dimension, 
rotational symmetry, and parity transformation properties, 

Um x ^ y < M*)MV) > = n r ^ 1 + 2r 2 {X ~ y)i{x ~ y)j ; (20) 

\x — y\ \x — y\ 6 

Um^ y - 2 < E t {x)E 3 {y) > = qxy^ + 2g 2 {x " v)i{x ~ y)j ; 
e 2 \x — y\ z \x — y\* 

$ij , nn {x-y)i(x-y)j 



lim x ^ y < Ei(x)Aj(y) > = s v . t= + 2s 2 - 

\x — yf \x — t/| 4 



en „ (x — y)i(x — v)i — (x — y)i(x — y)i 



\x — y\ \x — y 



The computational details are published elsewhere |T^] . Here, due to limitations of 
space, we only present the results, 

J = = -diEi (21) 
e 

J ee: := -2 : e ij A i E i : (22) 

The calculation of the spatial components of the current is more involved. The result 
that we get by expanding the defining relation eq.(|16p in powers of 1/A is, 

Ji = thfil> = -enAAi + \\: PE^AE) + jE 2 Ai :1 (23) 

eA L J 

, where E^ = e^Ej. The leading term is proportional to the UV cutoff. We therefore 
have to keep also the next order term, which vanishes in the naive continuum limit. 
(The cutoff-independent terms avearage to 0). The constants k, /3, and 7 depend on the 
regulator ratio |f/|/|e|. However, the final result should not depend on this ratio. This 
means, that the formally small terms that we discarded do not disappear completely, but 
renormalize the constants in eq . (123|) P1. The renormalized values of the coefficients can be 
determined by requiring that the currents satisfy the tree - level current algebra, possibly 
modified by Schwinger terms, and terms that vanish in the continuum limit, and by the 

3 One can show explicitly that the higher order terms lead to a finite renormalization of the coefficients 
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requirement of the restoration of Lorentz invariance in the continuum limit. The situation 
here is analogous to that in 1+1 dimensions [Q, where the use of Lorentz invariance is 
required to fix the overall scale of the current. The current is then determined as |L4], 



e/cA 



2E t (AE) - E 2 A/\ + —Et 
J e 



(24) 



The constants k and M will be given below. The fermionic bilinears are local functions 
of Ei and Ai. Their algebra can be calculated to give, 



[Ji{x), Jj(y)] = 2ie i:j J(x)5 2 (x - y); [Ji(x),J(y)] = -2ie ij J j (x)5 2 (x - y); 



(25) 



[J (x), Ji{y)} = -i 



Ak6 « + i^U 

2i 



(EiEj - 1/2E 2 6 K 



[J(x), J (y)} = -E i (x)e l3 d x J 5 I (x - y) 



The first two commutators are the canonical tree - level ones. The other two exhibit 
explicitly the Schwinger terms, which appear also at the one loop level in perturbation 
theory. 

We next calculate the bosonized form of the energy-momentum tensor of the theory. 
The gauge invariant, symmetric energy-momentum tensor of QED 3 is given by, 



rpflU rpfiV , Ji/lf 

13 ' F 

where the gauge and fermionic contributions are, respectively, 

1 



Tn" = F^Fx + -q^F 2 



and, 



I F 



(26) 



(27) 



(2f 



and where D = d — ieA is the covariant derivative. 

The general form of the bosonized Hamiltonian density following from eq. (|T^) is, 

T 00 = \b 2 + \e 2 + -|- (m) 2 + b^-AA 2 + \[ C :{A-E) 2 :+d: E 2 A 2 ]+fA-E (29) 

If the theory is to be Lorentz invariant with standard transformation properties for 
Ei,B,jQ, and Jj, Maxwell's equations should be satisfied in the continuum limit. We 
therefore require, 

B(x) = i[H, B{x)\ = -e^dtEj + 0(1/ A) (30) 

Et = i[H, E^ = -Ji + e^djB + 0(1/ A 2 ) (31) 

Since Jj itself contains terms of order 1/A, we demand that the second equation be 
satisfied to order 0(1/ A 2 ). 
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It turns out that this procedure still leaves the coefficients k and a undetermined. 
Our next step will be therefore to impose the complete Poincare algebra. 

A sufficient condition for Lorentz invariance of a theory invariant under spatial rota- 



tions and translations is the following commutation relation fl6| . 

-i[T 00 (x),T 00 (y)) = -(T?(x)+T?(yW6 2 (x-y) 



(32) 



It can be easily verified (by multiplying eq. (|32"D by linear functions of x and y and inte- 



grating) that eq. (|32|) implies the closure of the Poincare algebra []16 |, with the following 
genertors: the energy j d 2 xT 00 (x), the momentum Pj = / d 2 xT®, and the angular mo- 
mentum and boost generators, given, respectively by, 



^ = J d 2 x[x T° - XiT 00 }; L = J d 2 xe ij x i T° 



(33) 



We now require that the angular momentum generator has the standard "orbital'' 
and "spin" parts, 



L = d 2 xeijXiTf = / d 2 x[e ij x i E k d j A k - A ■ E] 



(34) 



Calculating the commutator in eq. (|3"2"D , and performing the operator product expan- 
sion using the forms of the propagators in eq.(pl| we find that Lorentz invariance is 
restored, in the limit A — > oo for, 



r 



AO 



2 B + 2 E + 2e 2 nA 
1 



im) 2 + *-kka 2 



(35) 



2kA 



-2(A ■ Ef + E 2 A 2 



+ MA-E: 



T = BE, - (djE^Ai 



2< .f > d i (AE) + ¥-d,(AiAi + AjAi) 
A s k 8 



(36) 



with, 



Ji 



-enkAi + 



en A 



2EjA ■ E - AjE z 



and k determined as, 



M = < A-E > 



A < E 2 > 



e 2 <B 2 > 



e 



(37) 
(38) 



We have now determined all the coefficients in the Hamiltonian density and the cur- 
rent in terms of ffg^ and < A-E >. We want to stress here that, although eq.(|3~5"D con- 
tains a parameter n, it does not define a one - parameter set of theories. The bosonized 
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version of QED3 corresponds to a unique choice of k, which satisfies eq.(^). Unfor- 
tunately, since eq.(|35D defines a strongly interacting theory, we cannot determine the 
numerical value of k. This would involve the solution of the model (at least in the UV 
region) for arbitrary k, calculating < B 2 > K and < E 2 > K , and solving the selfconsistency 
equation 



We have therefore expressed completely the one flavor QED 3 in terms of the electric 
fields Ei, and their conjugate momenta A\. The resulting theory is local, and Lorentz 
invariant. The fermionic bilinears satisfy a current algebra, which is modified from its 
tree level form by Schwinger terms. 

We now make several final comments. First, the fermionic operators eq.(|T3|) anti- 
commute only at distances larger than the ultraviolet cutoff. In fact, the phase in the 
anticommutation relations contains a power tail of the form (l/(A|x — y\) s , where s is a 
number of order one. Therefore, also in the lattice version of this theory, fermionic op- 
erators do not have canonical local anticommutation relations. This is the way in which 
the lattice theory avoids the standard doubling problem in the continuum limit. This 
aspect is different from the 1+1 dimensional case, where a theory of one staggered lattice 
fermion leads in the continuum limit to a theory of one Dirac fermion. Accordingly, the 
anticommutation relations in 1+1 dimensions are canonical at all distances up to the 
cutoff. 

Second, the fermionic operators we have constructed solve the Gauss's constraint. If 
one relaxes this condition, there are additional possibilities. One natural modification is 
to substitute for Ei in eq.flTTD the linear combination, H = Ei + Ce 2 6ijAj. This would 
give i[)>i{) = -diEi + CeB. However, the coefficient C in this case is not arbitrary. It 
is crucial for the consistency of the derivation that / d 2 xdjli has quantized eigenvalues 
TM. Since the magnetic flux in QED3 is quantized in units of — , the coefficient k must 



be given by ^- with integer n. This possible modification is a reflection of the well known 



regularization ambiguity in fermionic QED 3 ||17|| , which leads to the appearance of an 
induced Chern - Simons term in the action. 

Third, although the theory we have derived is Lorentz invariant, it is not written in 
terms of Lorentz covariant fields. The appearance of the UV cutoff in the expressions 
for the current and the Hamiltonian also point to the fact that the variables Ei and Ai 
are not the most convenient ones. It would be desirable to make one more step, and find 
a formulation in terms of a different, covariant bosonic field. We believe, that the most 
natural candidate for this field is the magnetic vortex field. In 2+1 dimensions it is a 



scalar fll3| . When written in terms of this field, the electric current is trivially conserved, 
and the electric charge has the interpretation of a topological winding number (which is 
quantized classically) [ p~8|j . 

Finally, we note that the fermionic operators which we have constructed involve a 
product of vortex and antivortex operators, carrying opposite magnetic fluxes. Therefore 
the operators do not carry net magnetic flux, and thus are not analogous to 2 + 1 
- dimensional dyons. The mechanism of their anticommutativity is different from that 
discussed in [[19[], and subsequently extensively exploited in the analysis of Chern - Simons 
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theories. In fact, our construction has a natural generalization to 3 + 1 dimensions, where 
the role of the vortex - antivortex pair is taken over by an infinitesimal closed t'Hooft 
loop of half - integer strength. Work along these lines is currently in progress. 
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